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Progress in modeling hypersonic 
turbulent boundary layers 

By O. Zeman 

1. Motivations and objectives 

A good knowledge of the turbulence structure, wall heat transfer, and friction 
in turbulent boundary layers (TBL) at high speeds is required for the design of 
hypersonic airbreathing airplanes and reentry space vehicles. This work reports 
on recent progress in modeling of high speed TBL flows. The specific research 
goal described here is the development of a second order closure model for zero 
pressure gradient TBL’s for the range of Mach numbers up to hypersonic speeds 
with arbitrary wall cooling requirement. 

2. Accomplishments 

In this report, new compressible models and theories that lead to their devel- 
opment are reviewed with the focus on compressibility effects in quasi-equilibrium 
turbulent boundary layers. The primary purposes are to report on a new second 
order closure model (SOC) developed for hypersonic TBL’s and to present com- 
parison of model results with experiments in zero pressure gradient TBL’s up to 
freestream Mach number M e = 10.3. The following section is a modified and ab- 
breviated version of the paper of Zeman (1993). The model described in subsection 
2.2.2 is a new contribution. 

2.1 Introduction 

Recent renewed interest in high speed aerodynamics has led to new developments 
in theory, simulation, and modeling of compressible turbulence. Availability of 
direct numerical simulations (DNS) of basic homogeneous compressible flows has 
greatly facilitated the development of new models for compressible turbulent flows. 

In view of the recent DNS results and experiments, it now appears that in many 
flows of practical interest, the turbulence cannot be treated by the so-called anelas- 
tic models, where the variation of averaged density and pressure are accounted for 
but not their fluctuating fields. In the past three years, this realization has lead 
to development of a variety of new models which account for the effect of fluctuat- 
ing divergence (dilatation) on turbulence. In this paper, we shall focus mainly on 
the compressibility effects pertaining to TBL’s with zero pressure gradients (ZPG). 

The paper is organized as follows: in the following section, we present the back- 
ground and review of the current representation of the dilatational terms in the 
modeling equations. The subsequent sections highlight the new SOC model for su- 
per/hypersonic TBL’s, make comparisons with experiments, and conclude with a 
discussion of compressibility effects in high speed TBL’s and their consequences for 
the turbulence models. 
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2.2 Background, review 

As shown first in Zeman 1 , the Favre- averaged energy governing equations for 
homogeneous compressible turbulence, in the absence of any forcing, can be written 
as 2 

l-A r = -{e s + e d -U d ) ( 1 ) 

^- = (e s + ei - Ud)c~ l , (2) 

where q 2 /2 = uju j/2 is the turbulence kinetic energy and c p T is the mean enthalpy 
(Favre and Reynolds averages are denoted by tilde and overbar, respectively ) and 
6 8 = vUjUJj is the solenoidal dissipation associated solely with the enstrophy UjUj 
(of the solenoidal velocity field). The compressibility effects are contained in two 
terms labeled e d and IL*. These are associated with the dilatational (or compressive) 
velocity field which has nonzero dilatation u h j (denoted hereafter by 6). Thus, the 
dilatation dissipation (,t = ^ i'0 2 and IT,; = p9 / p is the pressure-dilatation correlation 
(per unit mass). The compressive and solenoidal fields are strictly separable only in 
homogeneous turbulent flows; in TBL’s, the treatment of e d and e s as two distinct 
contributions to total dissipation is valid only approximately. 

2.2.1 Dilatation dissipation 

The need for a representation in turbulence models of dilatational dissipation e d 
associated with fluctuating Mach number has been now recognized by many authors 
(see e.g. Viegas and Rubesin 1991; Wilcox 1991). Computational results indicate 
that in TBL’s over insulated walls for M e < 9, the maximum values of M t are below 
0.3, and hence e d due to shocklet dissipation is insignificant. However, in hypersonic 
TBL’s with increasing wall cooling, the sonic speed a near the wall decreases and the 
M, levels grow larger. The shocklet dissipation then assumes a controlling role: it 
maintains Mj below a certain threshold level, which according to the computations 
is always below M ( = 1; we find this aspect of the dissipation physically appealing. 
The basic expression for the shocklet dissipation given in Zeman (1990) is 

U oe ?jF{M* t J<) (3) 

where t is a suitably defined turbulence lengthscale and F{ A/ ( * , A ) is a function 
of the r.m.s. Mach number. M* is related to the principal r.m.s. Mach number 

M, = ql\J-,RT through M* = J^M t . The parameter K is the kurtosis of 
the fluctuating speed intended to characterize intermittency of a particulai 

turbulent flow. The computed curves F(M* , A ) vs M* for different I\ have been 
given in Fig. 2 of Zeman (1990). It is of note that the dependence of e d on the 
specific heat ratio 7 ( through M * ) improves the correlation of mixing layer growth 
rate with M c , when the layer streams are gases with different values of 7 (Viegas 
and Rubesin 1991). 
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In modeling (3D) turbulence, the quantities q 3 /£ and e 9 are considered inter- 
changeable; however, it should be emphasized that (3) is valid also in 2D turbulence 
(DNS only) where typically q 3 /£ » e 9 . We also point out that no near- wall correc- 
tion is necessary in the expression for since F approaches zero at a much faster 
rate than the turbulent Reynolds number Rt (defined hereafter as R t = q A / 9e $ u). 
As in high speed mixing layers, F(M*,K) for TBL’s is approximated by an expo- 
nential function 

F(M;,K)) = ^ = c rf (l-exp{-(^L^) 2} ). (4) 

The parameters M to , and <tm are functions of the kurtosis K to approximate 

the shape of the F-curves for a specified K (see Zeman 1993 for details). 

2.2.2 Pressure- dilatation correlation in ZPG TBL's 

In inhomogeneous flows, nontrivial contributions to the pressure-dilatation term 
arise from the interaction between the mean density gradient Vp and fluctuating 
pressure field. The derivation of the density-gradient contribution ( pQ) p to the 
pressure dilatation has been presented in Zeman (1991, 1993). The form of the 
model for flat plate TBL’s (ZPG) is 


(p$) p = Tq 2 u\{p 2 ) 2 )z (5) 

In the SOC model, the contribution (5) is indispensable for assuring a proper (Van 
Driest) scaling of mean and fluctuating velocities in the inertial sublayer as shown 
in Zeman (1991, 1993). However, in the presence of wall heat transfer (cooling), 
the model (5) has proved to be ineffective in enforcing the correct scaling, and it 
had to be modified (as discussed briefly at the end of the following section). 


2.3 Closure of the compressible TBL equations 


In the boundary- layer approximation, the principal equations governing the mean 
flow field are the mass, momentum, and enthalpy conservation Favre- averaged equa- 
tions 



( 6 ) 


Pjjf = ~P,i ~ ( P u i u i ~ 2 P S *]),j (?) 

DT _ 

P~jj£ = -(pT'uj( 1 - F p ) - kTj)j + ~p{e s + e d - U d )c p 1 , (8) 

and the density is obtained from the equation of state p = p e /(RT + u\) where 
p € is the freestream (constant) pressure. F p in (8) denotes the ratio of pressure to 
enthalpy fluxes F p = pu x /{])c p T f Ui). In compressible TBL’s, F p is expected to be 
non zero, and Zeman (1993) proposed an expression 
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The coefficient (0.3) in (9) was chosen to recover the correct adiabatic wall temper- 
ature for the range 0 < M e < 11- The small M t limit, F p - M* is required by 
scaling arguments. 

It is of interest to note that if the turbulent fluctuations follow an adiabatic 
relation p oc then F p would be unity, and no heat would be transferred 

by turbulence? In the presence of heat sources, the compressive turbulence field 
is ineffective in transferring heat since it is virtually adiabatic; hence, the heat is 
transferred by the solenoidal turbulence only. In this sense, 1 — F p in (8) reflects 
the reduced mixing efficiency due to turbulence of acoustic origin. 

The remaining quantities needed to close the mean momentum and enthalpy 
equations (7) and (8) are the Reynolds stresses U{Uj and the heat fluxes T u,. 
General conservation equations for these quantities are shown in Zeman (1993). 
These equations contain the following terms requiring closure: pressure gradient- 
velocity and pressure gradient- temperature correlations denoted respectively bj II ij 
and II j, the triple-moment (transport) terms, the solenoidal dissipation e 9 , and the 
compressibility terms and IIj = pO/p. 

The pressure and transport terms are modeled in the same manner as their in- 
compressible (Reynolds averaged) counterparts and developed previously by Zeman 
and Jensen (1987) for atmospheric TBL’s (rough walls) and by Zeman (1990) for 
free compressible flows. Zeman (1993) modified the rapid part of II t; to account 
for the Reynolds number effects near smooth walls. This has been accomplished 
by making the coefficients associated with rapid terms, functions of the Reynolds 
number R t . For the asymptotically large values of R t , the rapid term coefficients 
converge to the values for the rough wall TBL as discussed above. 

The effect of the rapid pressure terms II, and II, , is best illustrated by writing 

closure equations for the shear stress U 1 U 2 and heat flux T l u 2 in 2D flat plate TBL. 
With £1 in streamwise and £2 in wall-normal direction, one obtains 


Du\U 2 _ ^ U\U 2 

Dt r 


0.4(t/ 2 — SAabntf 2 )fA,2 + T.T. 


( 10 ) 


£1^ = -C e — -(u 2 2 -5Aab n q 2 )f, 2 + T.T. (11) 

Dt r 

We can immediately see that apart from the transport terms (T.T.), (10) and 
(11) have a similar form which also suggests that |T'| oc |«i|. The similarity 
has been achieved by the novel formulation of the rapid part of II,. By neglect- 
ing the advection and transport terms, (10) and (11) reduce to algebraic relations 
u 7u 2 = — v T U\ f 2 •> and Vu 2 = -a r T a with the eddy viscosity and diffusivity v r 
and olt being proportional, i.e. 


vt oc ar oc tiz^[1 


— 5Ao 


&22 + 1/3 


( 12 ) 


Hence, in the algebraic approximation, the model yields a constant turbulent 
Prandtl number Pti — vt i the model constants were chosen so that Pv( 0.9, 
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the value which is supported by the DNS data in a channel flow, and by experiments 
in TBL’s. In ( 12 ), 6 ^ is the anisotropy tensor and Aa is a /Zrdependent coefficient 
in the rapid pressure model (described in Zeman 1993). 

The model equation for e g has a conventional form independent of M*, except 
for the wall treatment. The wall boundary value e 9 (y = 0 ) is determined from the 
approximate integral balance of the kinetic energy equation 

{P» - - Cd + n rf }dy = 0 

where mean convection is ignored and the no-slip condition has been used. Fur- 
thermore, to eliminate the unphysical wall singularity in the e,-equation and in the 
return-to-isotropy pressure terms (due to r = q 2 /e s — ► 0 ), the minimum r is set by 
the Kolmogorov time scale 

t - 5 \/f ’ <i3) 

as suggested by Durbin (1991). 

2.3.1 Modification of p9 in the presence of wall heat flux 

As mentioned earlier, in the presence of wall heat transfer, the model (p9) p in 
(5) has to be modified since the wall heat flux induces an entropic temperature 
field, giving spurious contributions to p6 p . Zeman (1993) proposed to decompose 
the temperature field on the adiabatic contributions T a (x,£) (corresponding to an 
insulated wall TBL) and on the entropic contributions T 3 (x,t) which arise due to 
the surface heat flux alone (no dynamic heating); the actual temperature field is 
the sum T = T s -f T a . The appropriate density gradient to be applied in ( 5 ) must 
be based on T a , i.e. p w —(T a ), 2 /T. The details of the determination of the 
entropic and adiabatic temperatures are presented in Zeman (1993). 

2.4 Comparison with boundary layer experiments 

The TBL computations are made by forward integration of the model equations 
starting with some initially thin TBL with the momentum thickness Reynolds num- 
ber Re$ = Ut9/v e = 200 — 500. The numerical scheme utilizes the compressible von 
Mises’ transformation (Liepmann & Roshko 1967) in the inertial and outer region of 
the TBL where y + > 30, and, in the region below y + = 30 (where advection terms 
are negligible), the TBL is solved as a parallel flow. The vertical velocity in this 
region is nonzero (due to density variation) and is eliminated by the transformation 
to 77 = Jq p/ p w dy. This computational method is effective and accurate. 

Fig. 1 is a sample of the model-experiment comparison: the streamwise r.m.s. fluc- 
tuations are plotted in a similarity form ( pu Vr w ) 1/2 vs. y/6 for a variety of Mach 
numbers and cooling rates. The universal behavior of the computed profiles is quite 
surprising; the cross-hatched area represents the measurements as compiled by Dus- 
sauge and Gaviglio (1987) (the first of these authors has pointed out to me that the 
low hot-wire value of u-fluctuations near the wall are likely due to errors associated 
with the transonic flow regime). 
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Figure 1 . Similarity profiles of the streamwise r.m.s. fluctuations; cross-hatched 
area represents the scatter of experimental data compiled by Dussauge and Gaviglio 
(1986). 

From a practical viewpoint, the most important test of a TBL model is its ca- 
pability to predict the friction coefficient Cf and the Stanton number St, defined 
as 

C/ " P'U? Cp Pe U e (T w - T aw y 

where q w is the wall heat flux and T aw , T w are the adiabatic recovery and actual 
wall temperatures. 

Fig. 2 shows standard plots of the ratio C f/C j 0 as a function of A/ e and T w /T aw 
where C fo is the low-speed value of C/(M e « 0) corresponding to the same Reg. 
Fig. 2a shows the model-computed values of Cf/Cf 0 vs M e , for an insulated-wall 
TBL for Re e » 10 4 and the Van-Driest II curve; a few data points are shown in the 
hypersonic range. Fig. 2b shows C f/C f 0 vs Tw/Taw for different M t . The unknown 
value of Cfo is assumed C fo = 0.02632I?e^° 25 . The model-computed values of Cj 
are in good agreement with theory and the data. 

The model-experiment comparisons of St vs C/ is shown in Fig. 3. The model 
values (for the range T w /T aw = 0.2 -0.6) indicate the Reynolds analogy factor Fr = 
2 St/Cj ~ 1.2; the displayed experimental values are in the range Fr = 0.9 — 1.2. In 
his review of experiments, Bradshaw (1977) suggests Fr be in the range 1.1 - 1.2. 
In view of the likely experimental errors, the model predictions of the principal 
parameters Cf and St are consistent with the data and theory. 

Fig. 4 consists of examples of the temperature profiles in the hypersonic range 
of Mach numbers. Fig. 4a shows model-experiment comparison of T/T e (y/6) for 
an insulated- wall TBL in helium with M e = 10.3; Fig. 4b is for M = 8.2 in air 
and with significant wall cooling (T w /T aw = 0.28). The computed temperature 




FIGURE 3. Variation of Stanton number St with C j for T w /T aw = 0.2 — 0.6. Ex- 
perimental data are from Laderman and Demetriades (1974), Kussoy and Horstman 
(1992), and Marvin and Coakley (1989). 
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profiles compare well with experiments (conducted on a sharp-edge flat plate). Of 
particular note is the prediction of the recovery temperature in the adiabatic TBL 
and of the temperature maximum near the wall in the cooled TBL. 

To demonstrate the model performance at low speeds, in Fig. 5 the modeled ve- 
locity profile (7 + (t/ + ) is compared with the DNS results of Spalart (1988). Although 
Spalart’s TBL Reynolds number, Re e = 1410, is below what is considered a min- 
imum self-similarity value, Re 0 = 3000, the model prediction is evidently in good 
agreement with the DNS. 

2.5 Model transition-to -turbulence prediction 

There are two kinds of transition to turbulence, what F. T. Smith calls the civ- 
ilized and the savage. In the civilized transition, small disturbances grow in ac- 
cordance with the appropriate instability mechanisms, eventually reaching a point 
where transition to turbulence is initiated by strong nonlinearities and formation 
of turbulent spots in the flow. In the savage, or by-pass transition, the stage of the 
orderly disturbance growth is bypassed, and turbulence is directly initiated by a 
nonlinear process. 

A fair indicator of the tendency to transition is the momentum Reynolds number 
Re$ of the pre-transition, laminar boundary layer. Ree accounts for the flow his- 
tory, and the transition Reynolds number Re$ t correlates well with the transition 
onset on flat plates. Typically, turbulence models use transition formulas which 
inform the model, on the basis of values of Re e , pressure gradient, and freestream 
turbulence intensity, when to turn on the eddy viscosity. Wilcox (1992) mentions 
the remarkable property of his k-u? model u to describe the nonlinear growth of flow 
instabilities from laminar flow into the turbulent flow regime.” In order to recover 
the appropriate transition Reynolds number for the Blasius profile, Wilcox modi- 
fied the model parameters (as functions of turbulent Reynolds number R*). Hence 
again, a correction has been provided to inform the model when to begin to amplify 
turbulence. 

A remarkable property of the present model is its capability to mimic transi- 
tion without any specific corrections added. This capability was tested only for 
high Mach numbers, and the results for two freestream Mach numbers are de- 
picted in Fig. 6. The computations started with a thin TBL with a relatively 
small Ree < 200, As seen in Fig. 6, the turbulence is initially attenuated and 
the TBL laminarizes. Only when Ree reached a certain (transition) value do the 
residual fluctuations within the boundary layer begin to rapidly grow until an equi- 
librium TBL is attained. More detailed investigations showed that (Ree)t increased 
with M t in a manner reminiscent of observed experimental transition (assuming 
(Ree)t oc y/Re x t ). 

It is of note that the transitional growth of turbulent energy first occurred in the 
upper part of the layer in the vicinity of the maximum of the mass vorticity pU iV 
(generalized inflection point). It is known from the stability theory that (pU iy )m a z 
is potentially a point of maximum instability growth. In the model, the coincidence 
between the maxima of q 2 and ]>U^y is a combined effect of the pressure- dilatation 
term p9 in (5) and of strong viscous damping near the wall. At Al e > 5, the 
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FIGURE 4 A& B. Model-experiment comparison of temperature profiles: a) data 
of Watson et. al. compiled in Fernholz and Finley (1977) under no. 73050504; 
M c = 10.31, R$ = 1.5 x 10 4 , in helium, b) data from Kussoy and Horstman (1992) 
with M e = 8.2, T w /T wa = 0.28, Re = 4.6 x 10 3 . 
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incipient maximum production of q 2 is always in the upper part of the layer near 
the generalized inflection point; after the transition, the point of maximum q 2 moves 
towards the wall. 

2.6 Discussion , Conclusions 

We have developed a new SOC model intended for general applications in high- 
speed turbulent flows. It incorporates the latest advances in compressible turbu- 
lence theories and modeling. The explicit compressibility effect on turbulence is 
represented by models for the dilatation dissipation and pressure dilatation term 
p6. Both €d and p6 depend on the r.m.s. Mach number ( M t ) and other structural 
parameters of the mean and fluctuating flow fields but not on the mean flow Mach 
number. The model predictions compare well with experiments for a wide range of 
Mach and Reynolds numbers. The importance of their contributions vary depending 
on flow speed and configuration. 

In some recently published work, the importance of explicit compressibility cor- 
rections in TBL models has been questioned. It is indeed possible to adjust incom- 
pressible models to perform well in compressible regime without compressibility 
terms. However, ours is a more fundamental question: are the compressibility ef- 
fects significant in reality, and can they be isolated in experiments and verified? If 
we consider DNS “experiments”, then the answer is obviously yes. Both the di- 
latation dissipation and pressure dilatation terms have been identified in DNS of 
shear-generated and rapidly compressed turbulence. Their seeming unimportance 
in TBL’s is only a question of degree. We find that as M e and wall cooling in- 
creases, ej becomes increasingly important. In the hypersonic regime with M e > 7 
and sufficiently strong wall cooling, the standard k-c models (without some form of 
dilatation dissipation) are likely to yield supersonic r.m.s. Mach number Mt > 1. 
This is obviously unrealistic; experimental evidence and DNS results suggest that 
Mt saturates well below unity. 

Concerning the importance of the pressure dilatation p6\ the density-gradient 
contribution to p6 constitutes a localized turbulence energy source which preserves 
the proper Van Driest scaling in the modeled TBL. The present results also suggest 
that p6 counteracts the damping, viscous effects which have a tendency to laminarize 
the boundary layer at high values of M t . The /^-contributions are also related to the 
ability of the SOC model to mimic transition to turbulence. We intend to address 
these matters and the plausibility of modeling transitional (high-speed) flows in 
future investigations. 

3. Future work 

We shall continue to refine the new SOC model and search out more data for 
model-experiment comparison. We also hope to apply the model in nonequilibrium 
situations such as a compression corner flow. 

In view of the ability of the SOC model to mimic transition, we shall attempt to 
investigate the connection between stability theory and model physics and to explore 
the potential of the SOC models to handle laminar and transitional regimes. 
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A major effort is going to be directed towards modeling nonequilibrium (incom- 
pressible) turbulence, such as in separated flows. 
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